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ABSTRACT 

Although the colour distribution of globular clusters in massive galaxies is well known 
to be bimodal, the spectroscopic metallicity distribution has been measured in only 
a few galaxies. After redefining the calcium triplet index-metallicity relation, we use 
our relation to derive the metallicity of 903 globular clusters in 11 early- type galaxies. 
This is the largest sample of spectroscopic globular cluster metallicities yet assembled. 
We compare these metallicities with those derived from Lick indices finding good 
agreement. In 6 of the 8 galaxies with sufficient numbers of high quality spectra we find 
bimodality in the spectroscopic metallicity distribution. Our results imply that most 
massive early-type galaxies have bimodal metallicity, as well as colour, distributions. 
This bimodality suggests that most massive galaxies early-type galaxies experienced 
two periods of star formation. 

Key words: galaxies: star clusters: general - galaxies: stellar content - globular 
clusters: general - galaxies: abundances 



1 INTRODUCTION 

Optical photometric studies of globular cluster (GC) sys- 
tems have shown that almost all massive galaxies have bi- 



modal GC colour distributions (e.g. Kundu & Whitmore 
2001 Larsen et al.p OOl). Since extensive spectroscopy has 
shown the majority of GCs are old (> 10 Gyr) (e.g. Forbes 
et al.|2001||Puzia et al.|2005||Strader et al.|2005[ ), this colour 
bimodality has usually been interpreted as a metallicity bi- 
modality. Metallicity bimodality suggests each galaxy has 



experienced two periods of intense star formation ( Brodie & 



|Strader||2006[ ) which must be accommodated into any sce- 
nario of galaxy formation. 

However, both |Richtler| ( [2006] ) and |Yoon et al."| (2006) 
showed that a strongly non-linear colour-metallicity relation 
can produce a bimodal colour distribution from a unimodal 
metallicity distribution. If the slope of the colour-metallicity 
relation flattens at intermediate colours, small differences in 
metallicity correspond to large differences in colour creat- 
ing a gap in the colour distribution where none exists in 
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the metallicity distribution. A unimodal GC metallicity dis- 
tribution therefore only requires one period of continuous 
GC formation. Thus the true form of the GC metallicity 
distribution therefore has important ramifications both for 
globular cluster formation and for galaxy formation. 

Although near-infrared (NIR) photometric studies (e.g. 
Kundu fc Zepf]|200"7| |Spitler et aT] |2008a| IChies-Santos] 



et al. 



2012) generally support the interpretation of colour 



bimodality as metallicity bimodality, large samples of spec- 
troscopic metallicities are required to settle the question of 
the form of GC metallicity distributions. Unfortunately for 
only a few galaxies have large numbers of GCs been studied 
spectroscopically. Although the Milky Way has been long 



known to host a bimodal metallicity distribution ( Harris & 



Canterna||l979 Zinn|1985 1, the GC metallicity distribution 
in M31 is less clear. |Barmby et aT| ( |2000| studied the colour 
and metallicity distribution of M31 GCs. They found colour 
bimodality in (V - K) but not in (V - I). Using 125 GCs 
with spectroscopic metallicities they found that the metal- 
licity distribution is bimodal. They showed that given the 
large reddening and photometric uncertainties, a unimodal 
colour distribution would be observed from their spectro- 
scopic metallicity distribution about 75% of the time. Using 
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Lick indices Galleti et al. (20091 derived the metallicities urates at metallicities greater than [Fe/H] ~ —0.5. Foster 



of 245 M31 GCs; they found that the metallicity distribu- 
tion is either bimodal or trimodal. Caldwell et al. (20111 



measured the metallicities of 282 GCs in M31 using the 
strength of iron spectral indices. Combining their spectro- 
scopic metallicities with 22 GC metallicities obtained using 
resolved colour magnitude diagrams, they found no statisti- 
cally significant metallicity bimodality (or trimodality) . 

|Beasley et al.| ( |2008[ ) measured the metallicity of 207 
GCs in NGC 5128, the closest, easily observable giant early- 
type galaxy, using Lick indices and found a bimodal metal- 
licity distribution. However the study of Woodley et al. 
( 2010bb, which measured the metallicity of 72 GCs in NGC 



5128 using Lick indices, did not find statistically signifi- 
cant metallicity bimodality despite finding bimodality in the 
metallicity sensitive [MgFe]' index and colour 



Alves-Brito 



et al. (20111 found evidence of metallicity bimodality in a 



sample of 112 GCs in the early-type spiral M104 (the Som- 



brero Galaxy) which they studied spectroscopically. Cohen 



et al. ( 1998) determined the metallicities of 150 GCs in M87 



using Lick indices and claimed evidence of bimodality. From 



this data set, Yoon et al. ( 2006J) noted that the metallicity 
sensitive index Mg 6 had a skewed single peak distribution 
rather than a bimodal distribution. Using 47 Lick index mea- 



surements by Cohen et al. (2003) of the massive elliptical 



M49's GCs, Strader et al. (20071 found metallicity bimodal- 



ity. Although these previous spectroscopic studies suggest 
GC metallicity bimodality is common in massive galaxies, 
a larger galaxy sample is required to confirm it. Due to the 
large amounts of telescope time required to acquire signif- 
icant samples of GCs, a more efficient observational tech- 
nique is required. 

The calcium triplet (CaT), at 8498 A, 8542 A and 8662 
A in the rest frame, is one of the strongest spectral fea- 
tures in the optical and NIR. In individual stars the CaT 
becomes stronger with both increasing metallicity and lower 
surface gravity (Ce narro et al.|2002 k Since giant stars dom- 
inate the NIR light of an old population and the surface 
gravity of metal rich giants is lower than metal poor gi- 
ants, the strength of the CaT in integrated light increases 
with metallicity. The CaT was first used by |Armandroff| 
& Zinn ( 1988 1 to measure the metallicities of Milky Way 



GCs using integrated light and has been successful at deter- 
mining the metallicities of individual stars (Rutledge et al 
1997 Battaglia et al. 2008). Single stellar population models 
( Vazdekis et al.|2003 hereafter V03) show that the CaT in- 
dex of integrated starlight is sensitive to metallicity while be- 
ing insensitive to ages greater than 3 Gyr. Although the CaT 
index is only weakly sensitive to the IMF for the |Salpeter] 
( 1955[ ) IMF and bottom light IMFs such as the |Kroupa| 
(20011 IMF (IV03I), it is affected by more bottom heavy IMFs 



such as those found in the most massive elliptical galaxies 
by |van Dokkum fc Conroy ( |2010[ ) and by |Cappellari et al.| 
( 2012[ ). 

Combining a wide field, red sensitive spectrograph (e.g. 
Faber et al.||2003 l and the CaT potentially allows 



DEIMOS 



large samples of spectroscopic GC metallicities to be assem- 
bled. Unfortunately the CaT has not yet been proven as a re- 



liable metallicity indicator. While Armandroff & Zinn ( 1988 1 



found a linear relationship between the strength of the CaT 
and metallicity over the range of metallicity of Milky Way 
GCs (-2 <[Fe/H]< 0), in the models of|V03]the CaT sat- 



et al. (2010 hereafter F10) was the first to use the CaT to 



determine the metallicity of a large number of extragalactic 
GCs. They used deimos to measure the metallicities of 144 
GCs in NGC 1407. Although they found a bimodal metallic- 
ity distribution, the fraction of GCs in the metal poor peak 
was different than the fraction of GCs in the blue colour 
peak. Their measurements of the CaT appeared to be non- 
linear with colour. They suggested this conflict is possibly 
due to saturation of the CaT, a non-linear colour-metallicity 
relationship or the effects of Paschen lines from hot hori- 
zontal branch stars on the CaT measurements. In addition, 
they found the brightest blue and red GCs to have simi- 



lar CaT values. More recently Foster et al. (2011 hereafter 



Fll) used the same technique to study the metallicites of 
57 GCs around NGC 4494. Unlike in NGC 1407, NGC 4494 
showed a linear relation between colour and the CaT in- 
dex. Although showing colour bimodality, the CaT strength 
appeared single peaked in NGC 4494. A larger sample of 
CaT measurements in multiple galaxies is thus required to 
establish whether the CaT can be reliably used to measure 
metallicity of extragalactic GCs. 



2 SAMPLE 

As part of the ongoing SAGES Legacy Unifying Globulars 
and Galaxies Survejj^] (SLUGGS, Brodie et al. in prep.), we 
studied the GC systems of 11 nearby (D < 30 Mpc) early- 
type galaxies. These galaxies cover a range of galaxy masses, 
morphologies and environments. Details of the galaxies in 
this study are given in Table [I] Further details of the obser- 
vations for NGC 3115 were presented in Arnold et al. ( 2011 1, 
for NGC 449 4 in|FlT]and for the remaining galaxies in |Pota| 
et al. (20121. We provide summaries of the data acquisition 



and reduction below. 



2.1 Photometry 



For most of the galaxies in our sample we used Subaru 



Suprime-Cam ( Miyazaki et al. 2002 1 gri imaging, supple- 



mented with other wide field imaging and HST imaging. 
The sources of our photometry is given in Table [2] We note 
that the g band imaging in NGC 821 and NGC 5846 and B 
band imaging from NGC 4278 were observed in poor seeing. 
All photometry is corrected for the effects of Galactic ex- 



tinction using the extinction maps of |Schlegel et all (|1998[) 



For galaxies and GCs that lack gri imaging we used em- 
pirical colour conversions to give all GCs equivalent (g — i) 
colours. These conversions are given in Appendix [A] Un- 
like our spectroscopy, the photometry used in this paper is 
heterogeneous - it was taken from different studies with dif- 
ferent zero points and aperture corrections. From the wide- 
field, ground-based imaging we selected GC candidates using 
cuts in colour-colour space. From HST imaging we selected 
GC candidates using cuts in colour-size space. We refer to 
the references for each photometric study for details of the 
selection. Colour magnitude diagrams and colour histograms 
of the photometric candidates are shown in Figure [I] 
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Table 1. Galaxy Properties 



Galaxy Hubble V sys (M — m) D M K A B N (g - i) 

Type (km s — 1 ) (mag) (Mpc) (mag) (mag) split 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



NGC 821 


E6 


1718 


31.85 


23.4 


-24.0 


0.48 


17 


0.97 


a 


NGC 1400 


E0 


558 


32.14 


26.8 


-23.9 


0.28 


34 


0.95 


b 


NGC 1407 


E0 


1779 


32.14 


26.8 


-25.3 


0.30 


202 


0.98 


b 


NGC 2768 


SO 


1353 


31.69 


21.8 


-24.7 


0.20 


19 


0.95 


b 


NGC 3115 


SO 


633 


29.87 


9.4 


-24.0 


0.21 


122 


0.93 


c 


NGC 3377 


E6 


690 


30.19 


10.9 


-22.8 


0.15 


84 


0.86 


b 


NGC 4278 


El 


620 


30.97 


15.6 


-23.8 


0.13 


150 


0.96 


b 


NGC 4365 


E3 


1243 


31.84 


23.3 


-25.2 


0.09 


131 


0.91 


b 


NGC 4494 


E2 


1342 


31.10 


16.6 


-24.1 


0.09 


53 


0.99 


d 


NGC 5846 


E0 


1712 


31.92 


24.2 


-25.0 


0.24 


54 


0.96 


b 


NGC 7457 


SO 


844 


30.55 


12.9 


-22.4 


0.23 


7 




c 



Notes Column (1): Galaxy name. Column (2): Hubble type from the NED database. Column (3): Systemic velocity from |CappelIari | 
et al.| ( |2011| . Column (4) and column (5): distance modulus and distance f rom|Cappellari et al.|(2011[ l based on surface brightness 
fluctuations (SBF) of |Tonry et~aT1 1 |2001| ). For NGC 4365 the SBF distance of |Mei et al.| | |2007| l is used. We assume that NGC 1400 and 

NGC 1407 are part of the same group and use the mean of their distance moduli. Column (6): K band absolute magnitude from 
|Cappellari et aL] | |2011[ l calculated from the 2MASS ( |Skrutskie et al.|2006[ ) total K apparent magnitude and the distance modulus in 

Column (4), corrected for foreground extinction. Column(7): B band foreground Galactic extinction from |SchlegeI et al.| l |1998| . 
Column (8): Number of GCs with CaT measurements. Column (9): (g — i) colour split between subpopulations from literature. NGC 
4365 shows evidence for colour trimod ality; see|Blom et al.| |2012b[ l. NGC 7457 appears unimodal; see |Hargis et all] ( |20 1 1 1 ) . Colour split 
references: a: |Spitler et""aT1 ( |2008b| >, b: |Pota et al.| ( |2012^ c: this work, d: |Fll| e: |Hargis et~al~1 \201l\ . NGC 1400, N GC 1407 and NGC 
3115 are not part of ATLAS dlJ sample so their radial velocities were taken from NED while their distances are from |Tonry et ah] j200l| 
with the same -0.06 mag distance moduli offset. Their K band absolute magnitude is also calculated from the 2MASS total K apparent 

magnitude and the distance in Column (3). 



Table 2. Photometry 



Galaxy 
(1) 



Instrument 
(2) 



Bands 
(3) 



Source 
(4) 



NGC 821 
NGC 1400 
NGC 1407 
NGC 2768 

NGC 3115 
NGC 3377 
NGC 4278 

NGC 4365 

NGC 4494 
NGC 5846 

NGC 7457 



Suprime-Cam 

Suprime-Cam 

Suprime-Cam 

Suprime-Cam 

HST ACS 

Suprime-Cam 

Suprime-Cam 

Suprime-Cam 

HST ACS 

Suprime-Cam 

HST ACS 

Suprime-Cam 

Suprime-Cam 

HST WFPC2 

WIYN 

Minimosaic 

HST WFPC2 



gri 
gri 
gri 
Riz 

BVI 
gri 
gri 

BVI 

gz 

gri 

9* 
gri 
gri 
VI 
BVR 

VI 



Pota et al. 


(2012} 


Spitler et al 
Spitler et al 


_ (12012} 

14113 


Pota et al.| 
Pota et al. 

- 


^20T2j 
IHIlL 


Arnold et al 
Pota et al. | 


J|201l| 


Pota et al. 


|20T2j 


Pota et al. 




Blom et al.| 


||20T2b| 


Blom et al.| 


\lWlb\ 


Fill 


Pota et al. 


J2012} 


Forbes et al 


rfT99ol 


Hargis et al.[(201lf 



Chomiuk ct al. ( 2008 ) 



Notes Column (1): Galaxy name. Column (2): Telescope and in- 
strument. Column (3): Photometric bands. Column (4): Reference 
to source of photometry. 



2.2 Spectroscopy 

Between 2006 and 2012, we used deimos on the Keck II 
telescope in multi-object mode to obtain spectra of large 
numbers of candidate GCs. These observations were pri- 
marily designed to obtain large numbers of GC radial ve- 
locities rather than to study GC stellar populations. Early 
observations used a central wavelength of 7500 A while later 



observations used a central wavelength of 7800 A. All obser- 
vations used the 1200 line mm -1 grating and 1 arcsec slits. 
This set up yields a resolution of AA ~ 1.5 A and covers 
the CaT region for almost all slits. We observed 2 to 11 slit 
masks per galaxy with typical exposure times of 2 hours per 
mask. 

The deimos data were reduced using the idl spec2d 
pipeline (Newman et al. 2012 Cooper et al. 20121. The 



pipeline performs flat fielding using internal flats, wave- 
length calibration using ArKrNeXe arcs and local sky sub- 
traction. In addition to the science spectrum, a fully prop- 
agated error array and a background sky spectrum are also 
produced by the pipeline for each object. After running the 
pipeline, heliocentric velocity corrections were applied. Ra- 
dial velocities were measured using the IRAF 2 \ procedure 
fxcor. To separate GCs from stars IPota et al. (20121) used 



the friendless algorithm of Merrett et al. ( 2003 1 for their 
galaxies. For NGC 3115 we used a velocity cut of 325 km -1 
to seperate stars and GCs while in the case of NGC 4494 
the GC system is well separated in velocity space from stars. 
We measured the signal-to-noise ratio (S/N) per Angstrom 
of the spectra using the mean S/N per pixel in the range 
8400 to 8500 A. 

To increase the number of GCs with CaT measure- 
ments that have literature metallicities, we supplemented 
our DEIMOS data with long slit spectra of two M31 GCs 
obtained using lris on Keck I ( Oke et al. 1995 I as poor 



2 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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weather targets on 2010 August 14. These observations used 
the 600/10000 grating, a central wavelength of 8448 A and 
a 1.5 arcsec slit. B012-G064 was observed for 3 x 300s while 
B225-G280 was observed for 4 x 300s. These lris observa- 
tions were reduced in the standard manner using IRAF. 



2.3 Comparison sample 

To compare our CaT metallicities with previous work we as- 
sembled a sample of literature metallicities based on Lick in- 
dices ( Worthey et al.|1994 1. We used the metallicity studies 
of|Cenarro et al.| ( |2007[ > for 17 GCs in NGC 1407, |Kuntschner 
et al.|(|2002[)"andNo rris (priv. comm.) for 14 GCs NGC 3115, 
Brodie et al.l (p005l) for 21 GCs in NGC 4365 andlChomiuk 



et al.| ([2008]) for 13 GCs in NGC 7457. All these studies de- 
rived metallicities in the same manner using % 2 minimiza- 
tion of Lick indices with the single stellar population (SSP) 



models of |Thomas et al.| (|2003[> a nd |Thomas et al.| (2004) 
(both hereafter |TMB03+TMK04l ) . We removed object 25a 
from Kuntschner et al. ( 2002 I as its colour is much bluer 



than the other GCs and its radial velocity is low enough to 
be a galactic star. For these four galaxies we used our own 
GC photometry. 

For M31 we used the metallicities from Caldwell et al.l 



(20111 who used an empirical relationship to derive the 



metallicities from the strength of Fe Lick indices. Since the 



Caldwell et al. (20111 metallicities are on the Carretta & 
Gratton ( 1997 1 scale, Eq uation [2] was used to convert the 
metallicities to the same TMB03+TMK04 scale used else- 
where in this paper. We used the M31 SSDS GC photome- 



try of Peacock et al. (20101. Due to the variable foreground 
and internal reddening of M31 we only used the 25 GCs 
with literature extinction values derived from resolved colour 



et al. (2005 


), Galleti et al. (2006), Mackey et al. (20071, Pe- 


rina et al. i 


20091, Caldwell et al. (20111, and Perina et al. 



(20111. For GCs in common with multiple studies we used 
the average value; however values from other studies were 



preferred to values from Caldwell et al. (2011). Rich et al 



(20051 provides two sets of reddenings: the values adopted 
and the values fitted as a free parameter. Except for B045- 
G108 and B311-G033 we used the adopted values; for these 
two clusters the fitted values provide better agreement with 
the colour-metallicity trend seen in other clusters. To ac- 
count for the large uncertainties associated with the red- 
dening correction we have added a 0.05 dex in quadrature 
to the photometric errors on the colours. In summary we 
have Lick index-based metallicities for 90 GCs in five galax- 
ies, of which we have measured the CaT strength of 31 GCs. 



3 CALCIUM TRIPLET MEASUREMENT 

The CaT lies in a region of strong sky lines. Although the 
DEIMOS pipeline does a good job subtracting the sky, the 
sky line residuals prevent line indices from being accurately 
measured directly on the observed spectra. To remove the 
sky line residuals we used the technique of |F10| and fit stel- 
lar templates to the spectra. The stellar templates, which 
cover a range of metallicities and a range of spectral classes 
(11 giants and 2 dwarfs ranging from F to early M), were 



observed using DEIMOS from 2007 November 12th to 14th us- 
ing an instrument set up comparable to the one used for the 
GCs. The templates were fitted using the pPXF pixel fitting 
code of ICappellari & Emsellem ( 2004 1 which also fits radial 



velocity, velocity dispersion and continuum shape. We fitted 
the 8425 to 8850 A rest frame spectral range and masked 
areas with strong sky lines from the fit. 

Before measuring the strength of the CaT on the fit- 
ted spectra, the continuum was normalised. Unlike in |F10| 
the continuum was normalised without human intervention. 
Instead, the continuum was modelled using a linear com- 
bination of Chebyshev polynomials of order zero to eight. 
Regions of the fitted spectra around known spectral lines 
were masked. A continuum model was then fit using least 
squares. Pixels more than 0.4 % below or 1 % above the 
model were rejected and the model recomputed. This pro- 
cess was repeated until no further pixels were rejected. Af- 
ter the initial fit the pixels were only masked if they devi- 
ated from the model. The same index definition (8490. OA 
to 8506.0A, 8532. OA to 8552.0A and 8653. OA to 8671. OA) 
as |F10| was used to measure the strength of the CaT on the 
normalised spectra. 

In IF10I the errors in the measured calcium index were 
estimated using the S/N of the spectrum. Although the S/N 
of the spectra dominates the errors, both metallicity and ra- 
dial velocity do have an effect on the errors. At high metal- 
licities stronger weak metal lines make it harder to fit the 
continuum while at certain radial velocities such as 200 km 
s _1 CaT lines fall on sky lines. In this work the uncertainties 
of the index measurements were derived by using a Monte 
Carlo resampling technique. One hundred realisations of the 
spectra were created from the fitted spectra and the ob- 
served error array. The spectral fitting, continuum normali- 
sation and index measurement process was repeated for each 
resampling. The 68th percentiles of the resamplings were 
used as confidence intervals. This procedure gives asymmet- 
ric error bars. After measuring the strength of the CaT and 
calculating the associated uncertainty, each spectrum was 
inspected visually to check the results of the automated fit- 
ting and normalisation. Spectra with large velocity disper- 
sions, a > 100 km s~ , were rejected as likely being contam- 
inated by background galaxy light. Spectra with poor sky 
subtraction where also rejected. The software used to mea- 
sure the CaT index and calculate errors is available from the 
authors upon request. 



4 ANALYSIS 

4.1 Metallicity scales 

Studies of extragalactic GC metallicities have used different 
metallicity scales. Since they include the effects of varying a 
element enhancement, the SSP models of |TMB03+TMK04| 
have been commonly used in extragalactic GC studies to 
derive metallicities from Lick indices. However studies that 
have used empirical relations based on the properties of 
Milky Way GCs tend to use the metallicity values from the 



Harris ( 1996 1 catalogue. Although originally based on the 



Zinn & West ( 1984 1 metallicity scale, the Harris catalogue 
has evolved with improved abundances and is now based 
( Harris|2010 1 on the Carretta et al. (20091 metallicity scale 
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(g-i) [Z/H]CaT (g-i) [Z/H]CaT 



Figure 1. Colour magnitude and metallicity magnitude diagrams for each galaxy. For each galaxy the upper left plot is a colour 
magnitude diagram, the upper right a CaT metallicity magnitude diagram, the lower left a normalised colour histogram and the lower 
right a normalised metallicity histogram. Black filled circles and solid lines are GCs with CaT metallicities; red small points and dotted 
lines are GC candidates. In the metallicity magnitude diagrams the centres of the points are colour coded blue or red depending on 
the photometric colour split. Blue GCs are metal poor; red GCs are metal rich. Median colour and metallicity error bars are plotted to 
the left and right. The horizontal dashed lines show Mi = —8.4, the turnover magnitude from [Villegas et al.| | [2010[ l, and Mi = —11, 
comparable to lj Cen's absolute mag nitude (Afj = -11.25, [Harris|l996| which corresponds to a mass of ~ 1 X 10 e Mq. On the colour 
magnitude diagrams the vertical dashed line is the colour split between red and blue subpopulations from Table ^ On the metallicity 
magnitude diagrams the vertical dashed line is colour split transformed into metallicity using Equation [To] Note that the metallicities of 
blue GCs remain on the metal poor side of colour division and vice versa. 
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Figure 1. Continued 



which uses large numbers (~ 2000) of high resolution spec- 
tra of individual GC stars. A challenged faced both by SSP 
models and empirical relations is that the Milky Way lacks 
the metal rich GCs seen in massive early-type galaxies. 

Since |TMB03+TMK04| models have been most com- 



monly used in extragalactic globular stellar population stud- 
ies, we used their metallicity scale in this paper. To place 
metallicity measurements and models that used different 
metallicity scales on the same scale as the one defined 
by the |TMB03+TMK04| models we derived empirical re- 
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lations using Milky Way GCs. Mendel et al. (20071 mea- 



sured stellar population parameters of Milky Way GCs using 
|TMB03+TMK04| models. We used their total metallicities 



([Z/H]) and the iron abundances ([Fe/H]co9) given by Car- 
retta et al. (|2009 \ to derive a linear relationship between the 



TMB03+TMK04| metallicity scale and the |Carretta et al.| 



et al. ( 2009 I between their iron abundance scale and the 



( 2009 1 iron abundance scale 

[Z/H] = 0.829[Fe/H] CO9 - 0.083 . (1) 

This relation has a rms difference of 0.138 dex. We combined 
the previous relation and the relation provided by |Carretta 



Car- 



retta & Gratton (19971 iron abundance ([Fe/H]cG97) scale, 
to give a relation between the Carretta fc Gratton| ( |1997[ ) 



iron abundance scale and the |TMB03+TMK04| scale 

[Z/H] = 0.943[Fe/H] CG97 - 0.085 . (2) 

This relation has a rms difference of 0.151 dex. We also 
derived a relation between the |TMB03+TMK04] metallicity 
scale and the scale defined by the 2003 version of the Harris 
catalogue. Using the |Mendel et alT] ( |2007| measurements we 
found the following relation: 

[Z/H] = 0.939[Z/H] HO3 + 0.043 . (3) 
This relation has a rms difference of 0.135 dex. 



4.2 Biases 

Before deriving metallicities we checked whether our CaT 
measurements were biased by velocity dispersion, radial ve- 
locity or S/N. We used V03 SSP model spectra, with an 



age of 12.6 Gyr and a Kroupa 



(2001 1 IMF, to create spectra 



with arbitrary velocity dispersion, radial velocity and S/N at 
each of the model metallicities. The IV03I models were used 
because they well match the reso lution of the deimos spec- 
tr a (~ 1.5 A for both). The |V03| models use the isochroncs 



of Girardi et al 



( 2000 1 and the empirical stellar library of 



Cenarro et al. (20011. With increasing velocity dispersion 



a bias towards lower index values was seen for the highest 
metallicities. However, this effect only became important for 
velocity dispersions larger than those seen in GCs. A bias 
towards higher index values was seen at certain radial ve- 
locities such as 200 km s _1 . This is likely caused by the 
CaT lines falling on masked skyline regions. A strong bias 
to higher index values was seen at S/N < 20 A -1 with the 
effect being more severe for metal poor spectra. The bias 
appears to reverse for the most metal rich spectra. The dif- 
ference is seen in the fitted spectra so it is due to the pPXF 
fitting rather than the normalisation process. To correct for 
this effect the mean CaT value was measured on spectra 
generated from the model spectra for a range of S/N ratios 
and metallicities. By using the mean index values at large 
(~ 300 A -1 ) S/N as the true index values, the following 
relation was used to correct for the S/N bias: 



CaT = CaT te:ctra „ + a x (S/N)" 

CaTtextraw < 4.96 

CaTt ea; t rau , > 4.96 

-0.343 C&TtextraW + 3.20 C&Ttextrau, < 5.80 
0.031 C'2LTt ex traw H~ 1-03 G^^textraw > 5.80 



— 22.5 C&Ttextraw ~\~ 122 
-3.33 G&Ttextraw H~~ 27.0 



(4) 
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Figure 2. Correction for the S/N bias in the CaT index values. 
The vertical axis is the CaT index value, the horizontal axis is 
the measured S/N. The small points are the mean measured CaT 
values from spectra generated from | Vazdekis et al.| |2003[ l model 
spectra at different metallicities and S/N ratios. The solid lines 
are the CaT values measured at high S/N for each model metal- 
licity. The dotted lines are Equation[4] a fit to difference between 
the measured points and the high S/N values. The hollow circles 
show the measured values after using Equation|4]to correct them. 
The colours in this plot scale with the metallicities of the model 
spectra used with blue being metal poor and red being metal rich. 
The S/N bais appears to be caused by pPXF. 



The S/N bias and the applied correction can be seen in 
Figure [2] All further CaT values in this study are corrected 
for the S/N bias. We restrict our analysis to spectra with 
a S/N greater than 8 A -1 . Using a S/N cut of 10 or 12 
A -1 in our analysis gave identical results so we chose to 
use the lower S/N cut to increase our sample sizes. Colour 
magnitude diagrams and colour histograms of the GCs with 
CaT measurements are shown in Figure [l] 



4.3 Repeated Measurements 

To assess the quality of the CaT index measurements we 
compare the 68 repeated measurements in Figure [3] The 
reduced x 2 value of the repeated measurements is 1.41 sug- 
gesting that that Monte Carlo resampling process is under- 
estimating the errors. Adding a systematic error of 0.255 A 
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Figure 3. Repeated measurements of CaT index values. The 
dashed line is one-to-one. The rms difference of the repeated mea- 
surements is 0.84 A. The error bars include the 0.255 A systematic 
error. 




Figure 4. Relation between the CaT index of |Armandroff fc| 
|Zinn| ([l988 ) and our index definition. Both indices are measured 
on spectra from the 1 2.6 Gyr,|Kroupa ||2001[ l IMF, single stellar 
population models of |Vazdekis et al.| 1 2003^ 7 The points are the 
measured values and the solid line is a quadratic fit (Equation^. 



in quadrature lowers the reduced % 2 to 1.00 We included this 
systematic error in future analysis. For repeated objects the 
simple mean of the CaT indices was used for further analy- 



4.4 Calibration 

Two approaches were taken to calibrate the relation between 
the CaT index and metallicity. First we used the same ap- 
proach as|F10| using the CaT index strengths of |Armandroff| 
& Zinn ( 1988 hereafter AZ88) to perform an empirical cal- 




Figure 5. Relation between metallicity ([Z/H]) and the CaT in- 
dex definition of |Armandroff fc Zinn] ||1988|l for 18 M ilky Way 
GCs. The [Z/H] values are from |Carretta et al.| l |2009[| converted 
to the same metallicity scale used in the models of |Thomas et ah] 
( |2003| > using Equation^ The solid line is a linear fit (Equation^ 
which has a rms scatter of 0.11 dex. 



ibration. We derived a relation between our CaT index and 
the |AZ88| index definition (CaTAzss)- We measured both 



our index and the I AZ881 index on[V03 SSP models with an 
age of 12.6 Gyr and a IKroupal pOOll IMF. The following 



quadratic relation was fit to the measurements and is shown 
in Figure [4] 

CaT AZ 88 = (-0.068 ± 0.007)CaT 2 + (1.55 ± 0.08)CaT 

+ (-0.96 ±0.22). (5) 

We also rederived the relationship between metallicity and 
thelAZ88lind ex. Unlike lFlOl who used the 2003 edition of the 
|Harris| ( |l996| l catalogue for metallicities, we used the metal- 



licities provided in Appendix A of Carretta et al. ( 2009 1 



converted to the TMB03+TMK04| scale together with the 
CaT AZ 88 values provided in Table 6 of |AZ88| Like |FT0| we 
ignore highly reddened clusters. Fitting a straight line gives: 



[Z/H] = (0.54 ± 0.02)CaT AZ 88 + (-3.90 ± 0.14) . 



(6) 



This relation is plotted in Figure [5] and has a rms scatter of 
0.11 dex. Combining Equations [5]and [6] gives: 

[Z/H] = (-0.037 ± 0.004)CaT 2 + (0.84 ± 0.06)CaT 

+ (-4.42 ±0.20) (7) 

which is plotted in Figure [6] This Figure gives the relations 
between our CaT index and metallicity for both of our cal- 
ibration approaches. 

Our second approach to calibrate the relation between 
metallicity and CaT index was to use the[V03]SSP. We mea- 
sured our index values on the model spectra with an age of 



12.6 Gyr and a Kroupa (2001 1 IMF and fitted a linear rela- 



tion between the metallicity of the models and the measured 
index value giving a relation of: 



[Z/H] = (0.461 ± 0.013)CaT ± (-3.750 ± 0.080) 



(8) 



which is also plo tted in Figure [6] T he rms of this relation is 
0.048 dex. Unlike |Tl^03±TMK04l|V03l does not correct for 
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Figure 6. Relation between metallicity ([Z/H]) and our CaT in- 
dex. The points are the index values measured on the 12.6 Gyr, 
Kroupa| (|2001| l IMF, single stellar population models of |Vazdekis| 
et al.| ( |2003^ . The model metallicities are corrected for the lo- 
cal stellar a element enhancement pattern using Equation [9] The 
solid line is a linear fit to these points (Equation^. For compari- 
son the dashed line is the relation derived using the measurements 
of |Armandroff &: Zinn| ( fl988] > (Equation [5| . Although both cali- 
bration techniques agree at high and low metallicities, there is a 
~ 0.3 dex difference in [Z/H] at CaT = 5.5 A. 



the effects of varying [cv/Fe] in the local stellar neighbour- 
hood. Following Mendel et al. (20071 we correct the model 



metallicities for this effect using: 



[Z/H] = 



[Z/H], 



[Z/H], 



< -1.0 



[Z/H]- 



- 1.28 



1.28 



-1.0 < [Z/H] V03 < 0.28 



[Z/H]. 



0.28 0.28 < [Z/H], 



(9) 



To check our results and to choose the better calibra- 
tion we compared the metallicities we derived with literature 
spectroscopic measurements. We restricted our analysis to 
GCs with a S/N greater than 8 A -1 and literature ages 
greater than 7 Gyr. The comparison of the CaT metallici- 
ties derived using the |AZ88] empirical calibration with litera- 
ture metallicities is given in Figure [7] For most metallicities 
([Z/H] < —0.5) this calibration gives metallicities 0.3 dex 
higher than literature values. However at higher metallici- 
ties there is good agreement. 

The comparison of the CaT metallicities derived using 
the |V03| SSP models with literature metallicities is also given 
in Figure [7] Unlike the empirical calibration, the SSP-based 
calibration approach gives metallicities in good agreement 
with literature values across the whole range of observed 
metallicities. Since the reduced \ 2 value is significantly lower 
for the SSP calibration (0.59) than for the empirical calibra- 
tion (2.26) and gives better qualitative agreement, we choose 
to use the SSP calibration approach for the rest of the anal- 
ysis (except when comparing to the previous measurements 
of F10| and |F11 in Section 5.4 1. Metallicity magnitude dia- 
grams and metallicity histograms for each galaxy are plotted 
in Figure [I] Colours, magnitudes, CaT index measurements 



and CaT metallicities for all measured GCs are given in Ta- 
ble [3] which is available in its entirety online. 

4.5 Colour— metallicity relation 

To compare the colours of GCs with their CaT metallici- 
ties we constructed an empirical colour-metallicity relation 



using literature metallicities (Section 2.3 1 and our own pho- 



tometry. We fitted the following broken linear function (See 
Figure M: 

(7.46 ± 1.28) x(g-i) + (-7.09 ± 1.01) 
(g - < 0.77 



[Z/H] 



(3.49 ±0.12) x (g 



(10) 



i) + (-4.03 ±0.11) 
(g - > 0.77 



The RMS of this colour-metallicity relation is 0.17 dex. Ver- 
sions of this relation in other optical colours are available in 
Appendix [B] Combining colour-metallicity data from mul- 
tiple galaxies may introduce scatter due to variations in the 
photometric zero points and reddening corrections. We note 
that the break in the colour-metallicity relation occurs at 
roughly the same metallicity ([Z/H] ~ —1.4) as the break 
in the metallicity-iron index relation used by Caldwell et al. 
(20111. 



A comparison of the metallicity distribution derived 
from the CaT compared to the distribution derived from 
the colours is given by the cumulative histograms in Fig- 
ure [9] To quantify this comparison we ran the Kolmogorov- 
Smirnov (KS) test on the CaT metallicities and the colour- 
based metallicities for the same objects. A non parametric 
test, the KS test measures the probability that two samples 
are drawn from the same distribution. We take a probability 
of less than p = 0.05 as evidence that the CaT-based and 
colour-based metallicities are drawn from different distribu- 
tions and inconsistent with one another. The results of the 
KS tests are given in Figure [9] 

4.6 Testing bimodality 

To test for colour and metallicity bimodality we used the 



Gaussia n Mixt ure Modelling (GMM) code of |Muratov fc 
|Gnedin| ( |2010 1. An improved version of the widely used 
KMM code of Ashman et all (119941), GMM calculates both 



the best fitting unimodal and best fitting bimodal distri- 
butions to the data. It then runs a parametric bootstrap 
to determine whether the bimodal distribution is preferred. 



Muratov & Gnedin (20101 discuss detecting bimodality and 



note that looking for an improved fit by two Gaussians over 
one is not so much a test of bimodality as a test of non- 
Gaussianity. Therefore further requirements must be met to 
declare a population bimodal. |Muratov fc Gnedin| ( |2010[ ) 
suggest using the relative separation of the fitted peaks, D, 
and the kurtosis of the sample to decide whether a sample 
is bimodal. They require D > 2 and negative kurtosis in 
addition to a significant bootstrap test result for a popula- 
tion to be bimodal. We adopt the same requirements. In line 
with previous bimodality studies we consider a distribution 
bimodal if bimodality is preferred over unimodality with a 
probability of p > 0.9. 

In addition to running GMM on the colours and metal- 
licities (both colour-based and CaT-based) of the GCs with 
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Figure 7. Left CaT metallicities derived using |Armandroff fc Zinn| ( |1988} based calibration compared to literature Lick index-based 
metallicities. Different galaxies and studies are denoted by different colours and symbols. The dashed line is one-to-one and the vertical 
dotted lines correspond to the highest ([Z/H]= —0.02) and lowest ([Z/H]= —2.02) mctallicity Galactic GCs used to derive the relation 
between index strength and metallicity. For metallicities of [Z/H]< —0.5, the metallicities given by empirical calibration are about 0.3 
dex higher than the literature values. The rms scatter is 0.27 dex and the reduced x 2 value is 2.26. Bight CaT metallicities derived 
using single stellar population (SSP) model of Vazdckis ct al. (2003) compared to literature metallicities. Different galaxies and studies 
are denoted by different colours and symbols. The dashed line is one-to-one while the vertical dotted line corresponds to the highest 
metallicity SSP model ([Z/H]= —0.05) used to derive the relation between index strength and metallicity. The lowest metallicity SSP 
population model used ([Z/H]= —2.32) is too low to appear on this plot. The rms scatter is 0.20 dex and the reduced \ 2 value is 0.59. 





NGC 4278 


J 7 ^ g 




150 GCs 




p = 0.128 




NGC 7457 


P k 




7 GCs 




p = 0.883 




[Z/H] CaT 
[Z/H] (g - i) 



-2 



-1 
[Z/H] 



Figure 9. Cumulative metallicity histograms for each galaxy. The solid black lines are the metallicities derived using the CaT; the 
dashed red lines are the same GCs with metallicities derived using (g — i) colours. The probability given by the KS test that the two 
samples are drawn from the same distribution is given in the lower right corner of each plot along with the number of GCs with CaT 
measurements . 
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Figure 8. Globular cluster colour-metallicity relations. Note that 
different studies used different techniques to derive metallicities 
and colours (see text for details) although we attempt place them 
on the same scales. Top Different colours and symbols refer to 
different galaxy studies. For [Z/H] > —1.4 the colour appears 
to be linear while for [Z/H] < —1.4 the relation steepens. Mid- 
dle Empirical colour-metallicity relations. The solid line is our 
broken linear fit to the data while the other lines are the colour- 
|Faifer et al.| ( |2011| l, |Sinnott et al.| | |2010| 
However, the different relations agree 
within the uncertainties of each relation above (g — i) ~ 0.8. 
We have measured CaT metallicities for the circled points. Bot- 
tom Single stellar population (SSP) model colour-metallicity re- 
lations. The solid line is our broken linear fit to the data. The 
coloured lines are the SSP mo dels of|Vazdekis et a l. ( 2010 ) , Maras- 
tonll|2005[l,|Yoon et aL] | |2011[ |, [C"onroy et al.| ( |2009[ | and |Raimondo 
et al.| l |2005| . All models are ~ 13 Gyr old; see text for details of 
the other stellar population parameters used. 



metallicity relations of 
and |Peng et al"] ( |2006 



calcium triplet measurements, we ran GMM on the colours 
of the GC photometric candidates. The photometric can- 
didates were restricted to those in the absolute magnitude 
range —13.5 < Mi < —8.4, the fainter limit corresponding 



to the turnover magnitude (Villegas et al. 20101 which is 



the peak of the GC magnitude distribution. We also limited 
the galactocentric radial range to radii less than that of the 
most distant spectroscopically confirmed GC. In the case of 
NGC 4494 we followed IF11I and excluded candidates closer 
than 0.5 arcmin to the galaxy. We did not make any correc- 
tions for contamination; such corrections would be different 
in each galaxy. We only ran GMM when the sample size was 
greater than or equal to 40. Results of running GMM on the 
metallicity distributions are given in Table [3] and on the 
colours in Table [5] The GMM metallicity results are plotted 
in Figure [T0"| 



5 DISCUSSION 

5.1 The calcium triplet as a metallicity indicator 

Figure [7] shows that metallicities measured using the CaT 
calibrated with SSP models agree with metallicities mea- 
sured using traditional Lick analysis better than the empir- 
ical, |AZ88f based, calibration. The agreement is seen over 
the range of observed metallicities and over a range of both 
GC magnitudes and galaxy masses. Since the CaT tech- 
nique relies on the strength of a single spectral feature, using 
the CaT to determine metallicities assumes that the age- 
metallicity, metallicity-calcium abundance and metallicity- 
horizontal branch morphology relationships are the same 
for all GCs in all galaxies as well as the same IMF. The 
CaT shares this issue with using colours or indices such the 
[MgFe]' index ( Thomas et al.|[2003 1 to derive metallicities. 

Unlike EOS we do not see the CaT saturate with metal- 
licity, at least for metallicities up to solar. We speculate that 
this is due to the effect of weak metal lines, some of which lie 
in the continuum pass bands of the CaT indices of |AZ88| and 
|Cenarro et al.| ( [2001[ ), acting to lower the pseudo-continuum 
of the CaT at high metallicities. Since the strength of the 
CaT increases and the pseudo-continuum decreases with 
metallicity, these indices saturate at high metallicity. As our 
technique fits the continuum, avoiding regions with weak 
metal lines, our continuum is less effected and our index 
does not saturate. 



5.2 The colour-metallicity relation of globular 
clusters 

In Figure [8] we compared our empirical colour-metallicity 
relationship (Equation 10 1 to those from literature (mid- 



dle panel) and to a range of different SSP models (lower 
panel). Comparing different relations is complicated due to 
the different metallicity scales used by various authors. We 
compared the following empirical relations: 



• The Faif er et al.| ( |2011[ ) relation was based on 
GMOS photometry and spectroscopy of GCs in NGC 3379, 
NGC 3923 and NGC 4649. The metallicities were derived 
using x 2 minimisation of Lick indices and the models of 
|TMB03+TMK04| 
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Figure 10. Comparison of colour and CaT-based metallicities for each galaxy. On the vertical axis is the CaT-based metallicity; the 
horizontal axis is colour-based metallicity. The plots on top are the colour-based metallicity histograms; to the right are the CaT-based 
metallicity histograms. The bimodal GMM fits, when available, are plotted over the histograms. The dashed line is a one-to-one line. 



The Sinnott et al. (2010) relation used photorn 



etry of NGC 5128 obtained using the Y4KCam imager on 
the Yale/SMARTS 1.0 metre telescope and metallicities de- 
rived using the strength of the [MgFe]' index calibrated to 
Milky Way GCs tied to the Harris metallicity scale. This 
relation was converted to the TM B03+TMK04| metallicity 
scale using Equation [3] 



• The |Peng eTaL] ( [20061 ) relation used ACS and SDSS 
photometry of M49 and M87 GCs together with unpublished 



photometry of Milky Way GCs. Metallici ties for M49 ( |Cohen 
|et al.|2003[ ) a nd M87 (|Cohen et al.|1998| > were derived using 
the models of Worthey ( 1994 1 and a limited number of Lick 



indices while Milky Way metallicities came from the [Harris] 
(19961 catalogue. The Peng et al. relation was converted 
from (g — z) to {g — i) using Equation Al and from the Harris 
metallicity scale to our metallicity scale using Equation [3] 

Although they disagree below (g — i) ~ 0.8, all three re- 
lations agree with our empirical colour-metallicity relations 
within the uncertainties of each relation. 

We also compared our empirical colour-metallicity rela- 



(2005), Yoon et al. (20111 



tion with the SSP models of Vazdekis et al. (2010 1, Maraston 



( 2009 1 and Rai- 



Conroy et al 

mondo et al.| ( |2005| |. We choose to compare the models at 
the same age (~ 13 Gyr). The metallicity scale of the mod- 
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Figure 10. Continued 



els is closer to the metallicity scale of Carretta et al. ( 2009 1 



than the |TMB03+TMK04| used in this paper so we used 
Equation [Tl to convert the model metallicities to the scale 



of[TMB03+TMK04. Below we briefly describe each set of 
models. 



• Version 9.1 of the SSP models of Vaz deklsl 
|et al.| ( |2010| ). We choose models with a |Kroupa| ( |2001[ ) 

IMF and with an age of 12.6 Gyr. Since these models only 
provide predictions in Johnson-Cousins filters, we use the 
transformation of |Jordi et al.| ( |2006[ > to convert (V — I) to 
(g — i). For the most metal poor GCs these models predict 
colours that are too blue while for metal rich GCs they 
predict colours that are too red. 



• The Maraston (2005 



SSP models. We choose 
models with a Kroupa (20011 IMF and with an age of 13 



Gyr. The Maraston models give a choice of blue or red 
horizontal branches at low metallicities; we choose the blue 
models as they provide better agreement with the data. 
Except for the most metal poor GCs these models predict 
redder colours at a given metallicity than most of the 
observed GCs. 



• Version 1.0 of unpublished Yonsei Evolutionary 
Population Synthesis (YEPS) models. These models 
are described in Yoon et al. (2006) and Yoon et al. (20111. 



Created to study the integrated colours of GCs, these mod- 



els assume a Salpeter ( 1955 I IMF and provide both scaled 



14 C. Usher et al. 




2.5 -2.0 -1.5 -1.0 -0.5 
[Z/H] (g - i) 



1.00 2 4 6 8 




1.0 



0.5 



0.0 



% -0.5 
u 

X 

N -1.0 



-1.5 



-2.0 



-2.5 



NGC 7457 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.00 
[Z/H] (g - i) 



Figure 10. Continued 



solar and a element enhanced models. We used models with 
[a/Fe] = 0.3 and an age of 13 Gyr. Although these models 
predict realistic colours at intermediate metallicities, they 
predict colours that are too red at high and low metallicities. 

• Version 2.3 of the Flexible Stellar Population 



Synthesis (FSPS) models of Conroy et al. (2009) 



We used 12.6 Gyr models and a Kroupa (20011 IMF. The 



FSPS models show good agreement with observations at all 
metallicities except at the metal rich end, where they are 
too red. 



• The Teramo SPoT models of Raimon do et al.l 
(2005). These models were used in the Cantiello & 



Blakeslee (20071 and Blakeslee et al. (20101 theoretical stud- 



ies of GC bimodality. These models assume a Scalo ( 1998 \ 



IMF and we used their models with an age of 13 Gyr. We 
converted the (V — I) colours provided by these models to 
(g — i) using the transform of Jordi et al. (20061. These 



models only provide colours similar to observations at in- 
termediate metallicities; at low and high metallicities the 
predicted colours are too red. 

Since we compare models with different adopted IMFs, 
we used the |Vazdekis et aL| ( |2010j ) models to investigate the 
effects of varying the IMF on the colour-metallicity rela- 



tion. Moving from a |Kroupa| (|2001) IMF to a [Scalo] ([1998 1 
IMF caused the (g — i) colour to become bluer but by less 



than 0.01 mag. Moving from a Kroupa (2001) IMF to a 
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Salpeter (19551 IMF caused the (g — i) colour to become SSP models of|V03Jand the colours from the SSP models of 



redder by 0.03 mag with the shift being independent of 



metallicity. IMFs more bottom heavy than Salpeter ( 1955 1 
have a larger red ward shift. Similar results were seen with 
the FSPS models. We conclude that the colour-metallicity 
relation is almost identical for bottom light IMFs such as 



heavy IMFs such as the Salpeter 



the Kroupa (20011 and the Scalo (1998) IMFs while bottom 



(19551 IMF have redder 



colours. Therefore the SPoT models which use a Scalo ( 1998 1 



IMF is directly comparable to the other models which use a 



Kroupa (2001 1 IMF while the YEPS models are likely a few 



hundredths of a magnitude too red due to their choice of a 



Salpeterj ([1955] IMF. 

All the models predict redder colours at the metal rich 
end than are observed. It is unknown whether this is due to 
a failing of the colour-metallicity predictions of SSP models, 
due to an issue with the ITMB03+TMK04I models used to 
derive the metallicities, or a mixture of the two. Out of the 
five SSP models we compare, the FSPS models of |Conroy] 
et al. ( 2009 I appear to provide the closest match to obser- 



vations. However it is difficult to quantitatively show that 
one model agrees with observations better than the others 
due to similarities of the models and the limitations of our 
data. 

In five of our eleven galaxies (NGC 1400, NGC 1407, 
NGC 2768, NGC 3377 and NGC 4494) the KS test (Fig- 
ure [9| provides evidence that the CaT-based metallicities 
are not drawn from the same distribution as the colour- 
based metallicities while in a sixth (NGC 4278) the KS test 
says there is a low but not significant probability that the 
two metallicities are not drawn from the same distribution. 
In the case of NGC 2768 the poor agreement between colour- 
based and CaT-based metallicity is likely due to the poor 
photometry available for this galaxy. The inner regions of 
NGC 3377 also suffer from poor photometry; however when 
GCs with poor photometry are removed, the disagreement 
between colour and CaT-based metallicity remains. 

Although NGC 1407 shows good agreement at high 
metallicities it displays poor agreement at low metallicities 
with the colour-based metallicities being more metal rich 
than the CaT-based metallicities (Fig ures [9}:, |lO|:) . NGC 
1407's GCs possess a strong blue tilt (Harris et al. 20061, 



a colour-magnitude (mass-metallicity) relation where the 
colour (metallicity) peak of the blue (metal poor) subpopu- 
lation becomes redder (more metal rich) with brighter mag- 
nitudes (higher mass) ( |Strader et al.|[2006| |Mieske et al.| 
2006). The blue tilt is thought to have caused by GCs self 



enrichment as more massive GCs are able to hold on to more 
of the metals formed by the first stars to form in the GCs 



(e.g. Bailin fc Harris||2009 1. It is plausible that this self en- 
richment changes the abundance pattern (for example the 
helium abundance) in such a way that the colour-metallicity 
relation (or at least the colour-CaT relation) is changed. The 
GCs we have CaT measurements for are effected by NGC 
1407's blue tilt as can be seen in Figure [If;. NGC 4278 is 
possibly another example of this effect as it also possesses a 
blue tilt (Usher et al. in prep.) and, like NGC 1407, shows 
weak evidence that the metal poor GCs are redder than in 
other galaxies (Figures[9j?, 10?). 

Another possible explanation is that the metal poor 
GCs in these galaxies have a different IMF than metal rich 
GCs in these galaxies or all GCs in other galaxies. Using CaT 



Vazdekis et al. (20101 we found that going from a Kroupa 



( 2001 ) IMF to a Salpeter ( 1955 ) IMF increased the measured 
CaT metallicity by 0.02 dex but make the {g — i) colour 0.03 
mag redder which due to the steeper colour-metallicity rela- 
tion would increase the metallicity derived from metallicity 
by 0.25 dex while a bottom heavy IMF similar that found by 
van Dokkum & Conroyl (120101) in the most massive ellipticals 



(slope x = —3) increased the measured CaT metallicity by 
0.07 dex but would make the colour 0.12 mag redder which 
would increase the colour based metallicity by 0.9 dex. Thus 
the redder colours of metal poor GCs in NGC 1407 and NGC 
4278 could be explained by a more bottom heavy IMF for 
these GCs than that of other GCs. 

Both NGC 3377 and NGC 4494 show the opposite 
trend, with agreement at low metallicities and the colour- 
based metallicities giving lower values than the CaT-based 
ones at high metallicities (Figures]!)]/', i, 10 f, i). In both NGC 
3377 and NGC 4494 the relationship between colour and 
metallicity appears to be linear and to be steeper than in 
other galaxies. The bluer colours of metal rich GCs in these 
galaxies could be explained by some of the metal rich GCs 
being younger than in other galaxies. In FSPS models of 
Conroy et al. (2009), at solar metallicity the (g — i) colour 



difference between a 7.9 Gyr model and a 13.3 Gyr model is 
0.10. Since the star formation history of early-type galaxies 



is dependent on mass and environment (e.g. |Trager et al. 


2000| |Terlevich & Forbes|2002 


Caldwell et al. 2003; Nclan 


et al.|2005 Thomas et al.|2005 


Smith et al. 2012 


i it is nat- 



ural to expect GCs in galaxies of different masses and en- 
vironments to have different age-metallicity relations. Some 
evidence of this is seen in the Milky Way where GCs associ- 
ated with dwarf galaxies that have been accreted on to the 
Milky Way show a different, steeper age-metallicity relation 
than those thought to have formed within the Galaxy it- 
self ( |Forbes fc Bridges|2010| potter et al.|2011[ ). NGC 1400 
is different again with agreement at low metallicities while 
having more metal rich colour-based values than CaT values 
at high metallicities (Figu res ]9f ), 1 10f >) . 

As noted in Section |2.l| our photometry is hetero- 
geneous. Different photometric zero points could shift the 
colour-metallicity relation along the colour axis or change 
its slope. We assumed one extinction value for all GCs in 
each galaxy so any differential reddening or internal ex- 
tinction could affect the colours of different GCs with in a 
galaxy to differing degrees. It is unclear whether our hetero- 
geneous photometry is responsible for the different colour- 
metallicity relations seen in different galaxies in our sample. 
If the colour-metallicity relation does vary from galaxy to 
galaxy it could be caused by variations in the ages or other 
stellar population parameters such as the IMF of the GCs 
between galaxies. 

5.3 Globular cluster metallicity distributions 

Although spectroscopic studies tend to sample the colour 
distribution well they are inherently biased to brighter mag- 
nitudes. If the stellar populations of the brightest GCs are 
different from the majority, then using only the brightest 
GCs will give a biased result. The blue tilt causes the metal- 
licity distribution of the brightest GCs to be more metal 
rich and less bimodal than the metallicity distribution of 
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et al.|2012b I and shows different kinematics for the different 
GC colour subpopulations (Blom et al. 2012a I so we also 



Figure 11. CaT metallicity distribution for NGC 4494. The solid 
line is a histogram of the measured metallicities. The dotted line 
is the two peak GMM fit while the dashed line is the three peak 
GMM fit. GMM returns a higher probability that three Gaussians 
are preferred over one (p = 0.90) than two Gaussians (p = 0.82) 
are preferred over one. This hints that NGC 4494's metallicity 
distribution is trimodal rather than unimodal or bimodal. 



the fainter GCs. We also note that the larger relative er- 
rors and smaller numbers in spectroscopic studies can make 
identification of metallicity bimodality more difficult than 
finding colour bimodality in photometric studies. 

Using GMM (Section [OJ Table [I] Figure [To}, we de- 
tected CaT metallicity bimodality with a probability of 
greater than 0.95 in NGC 3115, NGC 4278 and NGC 4365 
and greater than 0.90 in NGC 1407, NGC 5846 and NGC 
3377. NGC 821, NGC 1400 and NGC 7457 do not have 
enough measured GCs to determine whether their distri- 
butions are bimodal or not. In NGC 2768 we do not see 
metallicity bimodality nor is the colour distribution of GCs 
with CaT measurements bimodal. Due to the poor quality 
of the photometry for this galaxy it is unclear whether we 
should expect metallicity bimodality. 

NGC 4494 presents an interesting case. Although it does 
not meet one of our criteria for bimodality (a probability 
that two Gaussians are preferred over one of at least 0.9), 
GMM still gives a high probability that a two Gaussians are 
preferred over one of 0.819. The colour distribution of GCs 
with CaT measurements is bimodal (Table [5| so we would 
expect metallicity bimodality. Since the metallicity distri- 
bution of NGC 4494 appears trimodal (Figure 1 1 1 [ > we used 
GMM to test where the metallicity distribution is trimodal. 
We found peaks at [Z/H] = -1.73 ± 0.19, -0.69 ± 0.14 and 
0.26±0.11, widths of 0.17±0.12, 0.38±0.12 and 0.19±0.07 
and 10.5 ± 5.8, 29.5 ± 7.7 and 13.0 ± 5.0 GCs in each peak. 
GMM gives a probability of 0.895 that a trimodal distri- 
bution is preferred over a unimodal distribution, a higher 
probability than for bimodal over unimodal. The peaks are 
well separated (D = 3.55 and D — 2.32) and the kurtosis 
is negative k = —1.08 supporting the existence of trimodal- 
ity. The bimodal and trimodal GMM fits are plotted with 
the NGC 4494 metallicity histogram in Figure [IT] Due to 
the low number (53) of GCs these results should be treated 
with caution. IF11I found evidence for three kinematic sub- 
populations in that although the bluest GCs and the reddest 
GCs in NGC 4494 show evidence for rotation, the interme- 
diate colour GCs do not rotate. Taken together these lines 
of evidence suggest, but do not prove, that NGC 4494 has 
three GC subpopulations. 



ran GMM with three peaks on this galaxy. Although GMM 
finds peaks at [Z/H] = —1.53, —0.70 and —0.20 which are 
located at the same places as the colour peaks found by lBlom] 



et al. (2012b I transformed into metallicity, GMM does not 



find three Gaussians to be a significant improvement over 
one (p = 0.671). Due to the smaller numbers and larger 
errors of the CaT metallicities compared to the colours, it 
is not surprising that trimodality was not found to be an 
improvement. 

In general the colour distributions of the GCs with CaT 
measurements are similar to the photometric GC candidate 
colour distribution in eight galaxies with sufficient numbers 
for GMM to be run. However the blue peak of the GCs 
with CaT measurements in NGC 1407 is redder than the 
candidate blue peak in this galaxy due to the blue tilt while 
the red peak of the GCs in NGC 3377 is redder than that 
of the photometric GC candidates. The absolute magnitude 
limit of the GCs with CaT measurements varies dramatically 
as can be seen in the colour magnitude diagrams in Figure [l] 
In NGC 821, NGC 1400, NGC 1407 and NGC 5846, among 
the most distant galaxies in this study, the majority of GCs 
with CaT measurements are brighter than Mj = — 10 while 
in NGC 3115 and NGC 3377, the closest galaxies, GCs as 
faint as Mi = —8 have CaT metallicities. 

In Figure[l]; where the colour magnitude and metallicity 
magnitude diagrams of NGC 1407 are shown, the colour and 
metallicity distributions of GCs brighter than Mi = — 11 ap- 
pear different from those of fainter GCs. We split the GCs 
with CaT measurements in NGC 1407 into two magnitude 
bins and ran GMM on each. The results are plotted in Fig- 
ure [12] For GCs brighter than Mi = —11 the distribution 
looks quite unimodal. Fainter than Mi — —11 the distri- 
bution looks bimodal. GMM confirmed this with the bright 
population showed a low probability that 2 Gaussians are 
preferred to 1 (0.137) while the faint population showed 
high probability that 2 Gaussians are preferred (0.996). If 
only the 100 brightest GCs in NGC 1407 had been stud- 
ied spectroscopically bimodality would not have likely been 
found. This may explain why the Cohen et al.| ( |1998[ ) study 
of M87's GCs did not produce convincing proof of metallic- 
ity bimodality. It is known that in several galaxies including 
NGC 1407 ( |Romanowsky et al.||2009|>, M87 (Strader et aT 



[20TT] > and NGC 5128 ( |Woodley et~ah|2010a| ) GCs brighter 
than Mi ~ — 11 have different kinematics than fainter GCs 
further supporting the idea that the brightest GCs form a 
different population. 

Out of the eight galaxies with significant numbers of 
GCs with CaT metallicity measurements, we see bimodality 
in six and evidence for trimodality in another. In the eighth 
galaxy we do not see multimodality. However due to the 
poor quality of the photometry in this galaxy it is unclear 
whether we would expect bimodality based on the colour 
distribution of the observed GCs. Metallicity bimodality (or 
multi-modality) appears to be a common feature in massive 
early-type galaxies. 



5.4 Comparison with previous calcium triplet 
measurements 



NGC 4365 shows evidence for colour trimodality (Blom |F10| used the CaT to measure the metallicity of 144 GCs in 
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Figure 12. The CaT metallicity distribution of NGC 1407 split 
by magnitude. The black (red) histogram shows the distribution 
of GCs brighter (fainter) than Mi = — 11 while the dotted lines 
give the GMM fits. For the bright GCs GMM gives a low proba- 
bility that two Gaussians are preferred over one (p = 0.014) while 
for the bright GCs GMM gives a high probability (p = 0.996). 



NGC 1407 while |F11| used the same technique to measure 
the metallicity of 54 GCs in NGC 4494. In this study we have 
added additional GCs to the NGC 1407 data set and remea- 
sured the CaT for all GCs in NGC 1407 and NGC 4494. Un- 
like the previous studies we normalised the continuum with 
out human intervention (Section [3]| . Rather than calculate 
the uncertainty on the CaT measurements using the S/N as 
was done in |F10| and |F1 1| we used Monte Carlo resampling 
to calculate the errors which results in larger uncertainties. 
Unlike the earlier studies we corrected for an apparent S/N 



bias in the CaT index measurement (Section 4.2 1. We used a 



CaT metallicity calibration based on the SSP models of |V03| 
rather than an empirical calibration based on the strength 



of the CaT in Milky Way GCs as measured by |V03| (Sec- 
tion 



4.41 



Despite the somewhat different analysis we get qualita- 
tively similar results. The qualitative distributions of CaT 
indices in [FTo] and |FTl"1 both agree with this work, with 
NGC 1407 showing a bimodal distribution and NGC 4494 
not showing clear bimodality. The shapes of the colour-CaT 
index relations we find are similar to those of the previous 
studies with NGC 1407 showing the relation steepening at 
low metallicities and NGC 4494 showing a linear relation. 

In Figure [13] we compare our measurements to those of 
|F10| and |Fll| To be consistent with the calibration technique 
used in these papers we use the AZ88 metallicity calibration 
(Equation 7 1, rather than the SSP model based calibration 
(Equation 81, for this comparison only. Additionally we use 
Equation Rffin place of equation 2 of |F10| to convert the CaT 
values of |F10| and |Fll| into metallicities so any comparison 
is not biased by choice of metallicity scale. Although both 
measurements for NGC 1407 are consistent, for NGC 4494 
we find there is a significant offset with our measurements 
having metallicities 0.2 dex lower on average than those of 

Em 

Unlike [F10| we do not see significant Paschen lines in any 
of our fitted NGC 1407 spectra. We note that the strength 
of the Paschen line index used increases with metallicity due 
to the presence of weak metal lines in the index passband 
so most of the higher Paschen index GCs tend to be high 
metallicity GCs. Brodie et al. (20121 used SSP models to 
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show that horizontal branch morphology has little effect on 
the strength of the CaT. Like F10 we find that the brightest 



-1.5 -1.0 -0.5 0.0 
This Work [Z/H] CaT AZ88 

Figure 13. Top Comparison of CaT metallicities in NGC 1407 
from [Poster et alT] ( |2010| | with this work. Bottom Comparison of 
CaT metallicities in NGC 4494 from |Foster et aI.|j2TjTl) with this 
work. Over plotted with solid lines are linear fits to the metallici- 
ties while the dashed line is one-to-one. Both studies give consis- 
tent metallicities for NGC 1407. While the slope is consistent with 
unity, the NGC 4494 metallicities are 0.2 dex lower than those of 
this work. For consistency we use the |Armandroff fc Zin"n] | |1988| )- 
based calibration (Equation [TJl rather than single stellar popula- 
tion model metallicity calibration (Equation|8|l used elsewhere in 
this work and convert the metallicities of |Foster et al. | ( |2010^ and 
|Foster et aL] ( |2011} to our metallicity scale. 



GCs in NGC 1407 have similar CaT index values. However 
we find that the colours of the brightest GCs when trans- 
formed into metallicity are also similar. 



6 CONCLUSIONS 

We have used the strength of the calcium triplet (CaT) to 
measure the metallicities of 903 globular clusters (GCs) in 
11 galaxies ranging from brightest group ellipticals to iso- 
lated lenticulars. This is the largest sample of spectroscopic 
GC metallicities ever assembled and the first to provide 
large numbers (> 50) of metallicities for multiple galaxies. 
We showed that the CaT is a viable method of measur- 
ing the metallicity of extragalactic GCs with our CaT de- 
rived metallicities agreeing with literature Lick index-based 
metallicities. 

Using literature metallicities we denned a new colour- 
metallicity relation (Equation |10[). Although agreement is 
seen between the metallicities predicted by this relation and 
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CaT-based metallicities in half of our galaxies, the remaining 
galaxies show some disagreement either at low metallicity or 
at high metallicity. It is unclear whether these differences are 
caused by our heterogeneous photometry or real differences 
in the colour-met allicity relation between galaxies. Differ- 
ences in the colour-metallicity relation could be driven by 
differences in the age or the IMF of GCs between galaxies. 

We found spectroscopic metallicity bimodality in six of 
the eight galaxies with more than forty GC measurements, 
while evidence of trimodality was seen in one of the remain- 
ing galaxies (NGC 4494) and we can not rule out bimodality 
in the other (NGC 2768). We also found that the brightest 
GCs in NGC 1407 posses a different, unimodal metallicity 
distribution than the fainter, bimodal GCs. We caution that 
studies that rely on the properties of the brightest GCs to 
infer the properties of the GC system may give mislead- 
ing results. These bimodality results are similar to what has 
been seen in previous spectroscopic studies of the GC metal- 
licity distributions of early-type galaxies where three galax- 
ies (NGC 5 128, |Beasley et al.|2008| NG C 4594, |Alves-Brito 



et al.|201l| and M49, |Strader et al.|2007) show evidenc e of bi- 



modality while the situation in a fourth galaxy (M87 [Cohen| 
et al. 1998 1 is inconclusive. When the previous studies are 



combined with our own results, nine of the twelve galaxies 
show evidence for bimodality while in the remaining three 
galaxies bimodality can not be ruled out. 

The spatial distributions and kinematics of GCs sys- 
tem also favour multiple, distinct subpopulations. Red GCs 



are more centrally concentrated than blue GCs (e.g. Geisler 
et al.|[l996l |Bassino et al.|[2006| |Harris||2009l jFaifer et a" 



2011 



Forbes et al. 



2012). In addition the colour subpop- 



ulations often show different kinematics (Cote et al 



Schuberth et al. 



[2010] |Lee et al.||2010| |Strader et al" 



2003 



2011 



Pota et al. 2012 1 with sharp changes in rotation occurring 



at the blue/red dividing line. NIR photometric studies (e.g. 



Kundu fc Zepf|2007||Spitler et al.|20"08a||Chies-Santos et al.| 



2012 1 also see evidence for bimodality in several galaxies 



although there are unresolved issues with the connection 



between metallicity and NIR photometry (Blakcslcc et al. 
[2012] ). 

Our observed GC metallicity distributions, together 
with other observations of GC systems, favour a picture 
where most massive galaxies have bimodal GC systems. 
Since GC metallicity bimodality appears to be common most 
early-type galaxies should have experienced two periods of 
intense star formation. 
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Table 3. Calcium Triplet Measurements and Metallicities 
Name 



(1) 



RA. 

[deg] 
(2) 



Dec. 

[deg] 
(3) 



(9-i) 
[mag] 

(4) 



[mag] 
(5) 



CaT 

[A] 

(6) 



[Z/H] 

[dex] 
(7) 



NGC4494_GC100 
NGC4494_GC101 
NGC4494_GC102 



187.811375 
187.855092 
187.873400 



25.779614 
25.783281 
25.783786 



0.93 ± 0.03 
0.98 ± 0.08 
1.04 ±0.03 



21.04 ±0.02 
20.81 ± 0.06 
21.84 ± 0.02 



7 44+0.54 
6.52+ 1 - 47 



9.46 



0.54 
+0.64 
0.94 



pi oQ+0.25 
u -°^— 0.29 
-+0.68 
-0.25 

60+ - 30 
U.OU_ 4 3 



-0.75^ 



The full version of this table is provided in a machine readable form in the online Supporting Infromation. Notes Column (1): Globular 
cluster (GC) IDs. The IDs are the same as in |Pota et al.| | |2012[ l except for the GCs in NGC 4494 where 'NGC4494_' has been appended 
to the IDs from |FlT] Column (2) and (3): Right ascension and declination in the J2000.0 epoch, respectively. Column (4): Adopted 
(g — i) colour. Column (5): Adopted i magnitude. Column (6): CaT index measurement corrected for S/N bias. Column (7): Total 

metallicity. 



Table 4. Bimodal Metallicity Distribution Results using GMM 



Galaxy 


Sample 


N 


f-^poor 


&poor 


ftrich 


&rich 


fpoor 


V 


D 


k 


Bi 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 1407 


(s - 


202 


-1.01 ± 0.08 


0.19 ±0.04 


-0.27 ± 0.10 


0.38 ±0.05 


0.34 ± 0.11 


0.999 


2.50 ±0.28 


-0.96 


Y 




CaT 


202 


-1.52 ±0.26 


0.23 ±0.12 


-0.53 ±0.36 


0.52 ± 0.11 


0.16 ±0.22 


0.946 


2.46 ±0.66 


-0.40 


Y 


NGC 2768 


ig - 


49 


-0.64 ±0.28 


0.34 ± 0.10 


0.48 ±0.48 


0.11 ± 0.14 


0.96 ± 0.36 


0.645 


4.40 ± 1.51 


0.37 


N 




CaT 


49 


-1.12 ±0.44 


0.64 ± 0.20 


0.07 ± 0.43 


0.40 ±0.17 


0.78 ±0.29 


0.141 


2.24 ±0.68 


-0.78 


N 


NGC 3115 


(s - 


122 


-1.33 ±0.04 


0.32 ± 0.04 


-0.23 ±0.03 


0.23 ±0.03 


0.49 ± 0.05 


0.999 


3.94 ±0.37 


-1.06 


Y 




CaT 


122 


-1.19 ±0.10 


0.50 ±0.09 


-0.07 ± 0.06 


0.32 ± 0.05 


0.57 ± 0.08 


0.998 


2.67 ±0.50 


-0.75 


Y 


NGC 3377 


is - 


84 


-1.24 ±0.41 


0.49 ± 0.12 


-0.45 ±0.24 


0.16 ±0.12 


0.80 ±0.27 


0.903 


2.13 ±0.72 


-0.22 


Y 




CaT 


84 


-1.29 ±0.29 


0.54 ± 0.15 


-0.16 ±0.15 


0.35 ± 0.12 


0.56 ±0.17 


0.910 


2.49 ±0.62 


-0.66 


Y 


NGC 4278 


(s - 


150 


-1.08 ±0.07 


0.34 ±0.05 


-0.22 ±0.08 


0.25 ±0.05 


0.66 ± 0.09 


0.995 


2.86 ±0.37 


-0.82 


Y 




CaT 


150 


-1.48 ±0.15 


0.42 ±0.09 


-0.48 ±0.05 


0.32 ±0.06 


0.41 ± 0.11 


0.999 


2.65 ± 0.53 


-0.47 


Y 


NGC 4365 


(s - 


131 


-1.26 ±0.07 


0.21 ±0.04 


-0.36 ±0.06 


0.33 ± 0.04 


0.30 ± 0.06 


0.999 


3.26 ±0.31 


-0.97 


Y 




CaT 


131 


-1.47 ± 0.23 


0.25 ±0.10 


-0.51 ± 0.12 


0.40 ±0.07 


0.18 ±0.16 


0.951 


2.88 ±0.53 


-0.58 


Y 


NGC 4494 


ig - 


53 


-0.99 ±0.18 


0.35 ± 0.08 


-0.31 ± 0.12 


0.13 ±0.05 


0.71 ± 0.16 


0.954 


2.55 ±0.67 


-0.58 


Y 




CaT 


53 


-0.90 ±0.38 


0.62 ±0.19 


0.26 ±0.32 


0.19 ±0.17 


0.80 ±0.27 


0.819 


2.51 ± 0.53 


-1.08 


N 


NGC 5846 


ig - 


54 


-2.79 ±0.91 


0.14 ±0.32 


-0.63 ±0.33 


0.55 ± 0.16 


0.02 ± 0.30 


0.841 


5.37 ±2.20 


1.19 


N 




CaT 


54 


-1.42 ±0.31 


0.35 ±0.14 


-0.41 ± 0.13 


0.34 ± 0.08 


0.24 ±0.18 


0.925 


2.93 ±0.99 


-0.43 


Y 



Notes Column (1): Galaxy name. Column (2): Colour or CaT-based metallicity. Column (3): Number of GCs with CaT measured. 
Column (4): Mean metallicity of the metal poor subpopulation. Column (5): Dispersion of the metal poor subpopulation. Column (6): 
Mean metallicity of the metal rich subpopulation. Column (7): Dispersion of the metal rich subpopulation. Column (8): Fraction of the 
clusters in the metal poor population. Column (9): p-value that a bimodal fit is preferred over a unimodal fit. Column (10): Separation 

of the GMM peaks normalised by their width. Column (11): The kurtosis of the sample. Column (12): Bimodality of the sample. A 
p-value greater than 0.9, a separation of D > 2 and a negative kurtosis are all required for a sample to be considered bimodal. We did 
not run GMM on NGC 821 (17 GCs), NGC 1400 (34 GCs) and NGC 7457 (7 GCs) due to the low number of GCs with CaT 

measurements in these galaxies. 
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Table 5. Bimodal Colour Distribution Results using GMM 



Galaxy 


Sample 


N 


f-^blue 


^ olue 






£ 

J blue 


V 


D 


k 


Bi 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7? 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 821 


All 


110 


0.81 ± 0.05 


0.04 ±0.03 


1.03 ± 0.10 


0.16 ± 0.04 


0.18 ± 0.26 


0.969 


1.87 ± 0.70 


-0.48 


N 


NGC 1400 


All 


1010 


0.80 ± 0.01 


0.06 ± 0.01 


1.09 ± 0.02 


0.10 ±0.01 


0.53 ±0.05 


0.990 


3.53 ± 0.37 


-1.16 


Y 


NGC 1407 


All 


3312 


0.82 ± 0.00 


0.06 ±0.00 


1.12 ±0.01 


0.11 ± 0.01 


0.48 ± 0.02 


0.990 


3.36 ± 0.12 


-1.13 


Y 




CaT 


202 


0.86 ± 0.02 


0.05 ± 0.01 


1.08 ± 0.03 


0.11 ± 0.01 


0.34 ± 0.11 


0.999 


2.50 ± 0.27 


-0.97 


Y 


NGC 2768 


All 


195 


0.70 ± 0.04 


0.05 ± 0.02 


1.00 ± 0.03 


0.16 ±0.01 


0.17 ±0.09 


0.999 


2.51 ± 0.23 


-0.78 


Y 




CaT 


49 


0.97 ± 0.08 


0.10 ±0.03 


1.29 ± 0.14 


0.03 ± 0.04 


0.96 ±0.36 


0.645 


4.40 ± 1.51 


0.37 


N 


NGC 3115 


All 


180 


0.77 ± 0.01 


0.09 ±0.01 


1.07 ± 0.01 


0.07 ± 0.01 


0.54 ± 0.04 


0.999 


3.88 ± 0.26 


-1.24 


Y 




CaT 


122 


0.79 ± 0.01 


0.06 ±0.01 


1.09 ± 0.01 


0.07 ± 0.01 


0.47 ± 0.05 


0.999 


4.74 ± 0.39 


-1.38 


Y 


NGC 3377 


All 


71 


0.75 ±0.02 


0.05 ± 0.02 


0.94 ± 0.03 


0.07 ±0.02 


0.52 ± 0.13 


0.997 


3.07 ± 0.54 


-1.02 


Y 




CaT 


84 


0.79 ± 0.02 


0.08 ± 0.01 


1.01 ± 0.03 


0.06 ± 0.01 


0.67 ± 0.10 


0.994 


3.32 ± 0.46 


-1.01 


Y 


NGC 4278 


All 


700 


0.82 ± 0.01 


0.08 ±0.01 


1.10 ±0.01 


0.10± 0.01 


0.63 ±0.03 


0.999 


3.28 ± 0.20 


-0.82 


Y 




CaT 


150 


0.80 ±0.03 


0.05 ± 0.02 


1.02 ±0.04 


0.11 ± 0.02 


0.40 ±0.13 


0.999 


2.56 ±0.49 


-1.02 


Y 


NGC 4365 


All 


2159 


0.79 ± 0.00 


0.06 ±0.00 


1.05 ±0.01 


0.12 ±0.00 


0.36 ± 0.02 


0.990 


2.82 ±0.10 


-1.01 


Y 




CaT 


131 


0.80 ±0.01 


0.04 ±0.01 


1.05 ±0.01 


0.10 ±0.01 


0.28 ±0.06 


0.999 


3.34 ± 0.31 


-1.08 


Y 


NGC 4494 


All 


127 


0.86 ± 0.02 


0.09 ±0.01 


1.07 ±0.03 


0.03 ± 0.01 


0.80 ±0.09 


0.999 


3.07 ± 0.32 


-0.97 


Y 




CaT 


53 


0.86 ±0.02 


0.08 ±0.01 


1.06 ±0.02 


0.04 ± 0.01 


0.65 ± 0.11 


0.986 


3.20 ± 0.53 


-1.14 


Y 


NGC 5846 


All 


894 


0.76 ± 0.02 


0.07 ±0.01 


1.01 ± 0.02 


0.14 ± 0.01 


0.34 ± 0.08 


0.999 


2.29 ± 0.16 


-0.89 


Y 




CaT 


54 


0.89 ±0.09 


0.13 ±0.05 


1.05 ±0.13 


0.14 ± 0.05 


0.50 ±0.30 


0.022 


1.18 ± 1.14 


-0.24 


N 


NGC 7457 


All 


46 


0.85 ±0.04 


0.17 ±0.05 


1.43 ±0.18 


0.04 ± 0.08 


0.98 ±0.22 


0.551 


4.70 ± 1.49 


0.44 


N 



Notes Column (1): Galaxy name. Column (2): Whether all GC candidates or just GCs with CaT measurements. Column (3): Number 
of GCs. Column (4): Mean colour of the blue subpopulation. Column (5): Dispersion of the blue subpopulation. Column (6): Mean 
colour of the red subpopulation. Column (7): Dispersion of the red subpopulation. Column (8): Fraction of the clusters in the blue 
population. Column (9): p-value that a bimodal fit is preferred over a unimodal fit. Column (10): Separation of the GMM peaks 
normalised by their width. Column (11): The kurtosis of the sample. Column (12): Bimodality of the sample. A p-value greater than 
0.9, a separation of D > 2 and a negative kurtosis arc all required for a sample to be considered bimodal. 
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APPENDIX A: EMPIRICAL COLOUR 
TRANSFORMATIONS 

For galaxies and GCs that lack gri imaging we used em- 
pirical colour conversions to give all GCs equivalent (g — i) 
colours and i magnitudes. To convert HST ACS gz photom- 
etry into Suprime-Cam gi photometry we derived the fol- 
lowing relations using 169 spectroscopically confirmed GCs 
with both ACS and Suprime-Cam photometry in NGC 4365: 

(g-i) = (0.735±0.009)x( ff -z) + (0.147±0.012) mag (Al) 
and: 

i = z+(0.554±0.032)x(g- 2 ) + (-0.542±0.034) mag . (A2) 



These relations were used to convert the ACS gz photometry 
in NGC 4278 and NGC 4365. We also used the same objects 
to derive a relation between (g — r) and (g — i): 

(g-i) = (1.639 ± 0.011) x (g - r) + (-0.127 ± 0.007) mag . 

(A3) 

To convert (V — I) colours into (g — i) we derived the 
following relation using 30 photometric GC candidates in 
NGC 5846 brighter than i = 23: 

(g-i) = (1.259±0.076)x(\/-/)-(0.411±0.069) mag . (A4) 

To convert I to i we used the same data to get: 

i = I+(0.169±0.091)x(V-I) + (0.143±0.099) mag . (A5) 

We used these equations to convert the WFPC2 VI pho- 
tometry in NGC 5846 and NGC 7457. 

To convert BI photometry to gi we used the ACS BI 
photometry of |Forbes et al] ( [2006} of NGC 1407 and 57 
spectroscopically confirmed GCs in NGC 1407 to find: 

(g-i) = (0.651±0.015)x(B-/)+(-0.187±0.025) mag (A6) 

and: 

i = /+(0.162±0.146)x(B-/) + (0.081±0.259) mag . (A7) 

We used these relations to convert the ACS BVI photome- 
try in NGC 2768. Using 30 spectroscopically confirmed GCs 
in NGC 2768 with both transformed ACS and Suprime-Cam 
photometry we derived a relation between (r — z) and (g — i): 

(g-i) = (0.913±0.227)x(r-z) + (0.408±0.154) mag . (A8) 

We used this relation to convert the colours of GCs with 
only Riz Suprime-Cam photometry in NGC 2768. 

To convert BVR photometry into a gi photometry we 
used the WIYN Minimosaic BVR photometry of NGC 821 



from Spitler et al. ( 2008b I together with 42 spectroscopically 



confirmed GCs to get the following: 

(g - i) = (1.457 ± 0.101) x (B — R) + (-0.878 ± 0.103) mag 

(A9) 

and: 

i = R + (-0.565 ± 0.086) x (B - R) + (0.385 ± 0.084) mag . 

(A10) 

We used this to convert the WIYN Minimosaic BVR pho- 
tometry of NGC 7457. 

For NGC 4278 |Pota et al.| pOL2| prefer ACS photom- 
etry to Suprime-Cam photometry. When ACS photometry 
is unavailable they convert the Suprime-Cam BVI photom- 
etry into gz. We used Equations |A1| and |A2| to convert gz 
into gi. 
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APPENDIX B: EMPIRICAL 
COLOUR-METALLICITY RELATIONS 

Using the preceding colour-colour relations we converted 
Equation [To] into other optical colours: 

' (5.47 ± 0.94) x (g- z) + (-5.95 ± 1.02) 

(g - z) < 0.1 



[Z/H] = 



(2.56 ± 0.09) x (g - z) + (-3.50 ± 0.11) 
(g -z)> 0.84 



(Bl) 



[Z/H] = 



' (12.23 ± 2.10) x (g - r) + (-8.04 ± 1.01) 

(g-r) < 0.55 

(5.72 ± 0.20) X (g - r) + (-4.47 ± 0.11) 
(g-r) > 0.55 



' (9.36 ±1.71) X (V-I) + (-0.16 ±1.24) 

(V-I) <0.94 



[Z/H] = 



(4.39 ± 0.30) x (V - I) + (-5.46 ± 0.27) 
(V - 7) > 0.94 



(B2) 



(B3) 



' (4.86 ± 0.84) x (B — I) + (-8.48 ± 1.04) 

(B-I) < 1.47 



[Z/H] = I 



(2.27 ± 0.09) x (B - I) + (-4.68 ± 0.14) 
(B-I)> 1.47 



(B4) 



